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In this paper, we study topological phase transition in a wire medium operating at infrared 
frequencies. This transition occurs in the reciprocal space between the indefinite (open-surface) 
regime of the metamaterial to its dielectric (closed-surface) regime. Due to the spatial dispersion 
inherent to wire medium, a hybrid regime turns out to be possible at the transition frequency. Both 
such surfaces exist at the same frequency and touch one another. At this frequency, all values of the 
axial wavevector correspond to propagating spatial harmonics. The implication of this regime is the 
overwhelming radiation enhancement. We numerically investigated the gain in radiated power for 
a sub-wavelength dipole source submerged into such the medium. In contrast to all previous works, 
this gain (called the Purcell factor) turns out to be higher for an axial dipole than for a transversal 
one. 


I. INTRODUCTION 

Electromagnetic metamaterials are artificial media 
structured on the sub-wavelength level which exhibit ex¬ 
traordinary electromagnetic features and functionalities 
that are not found in nature. Most popular examples of 
these unusual functionalities can be found in works 13 . 
Many metamaterials possess unusual dispersion proper¬ 
ties. For regular metamaterials (lattices) their dispersion 
properties are described in terms of dispersion surfaces, 
also called Fresnel’s wave or iso-frequency surfaces. A 
dispersion surface defines wavevectors for all propagation 
directions at the given frequency in the reciprocal space 
of the lattice. One of the striking examples of metamate¬ 
rials are hyperbolic ones called this way because their dis¬ 
persion surface is hyperboloid. Such metamaterials are 
exploited in a variety of applications ranging jrom sub¬ 
wavelength imaging to radiative heat transfer®. Indeed, 
a non-magnetic hyperbolic metamaterial is most popular 
in optics, where magnetic properties are very difficult to 
realize. Optical hyperbolic metamaterials are anisotropic 
media whose principal components (more exactly - real 
parts) of the effective permittivity tensor have different 
signs. Because of the indefinite sign of the permittiv¬ 
ity tensor trace such media are also called indefinite 5 . 
Most expanded are indefinite (hyperbolic) media of uni¬ 
axial type. The different signs of the axial and transverse 
components of the permittivity tensor imply the hyper¬ 
boloid of revolution in the reciprocal space. Since the 
hyperboloid is infinitely extended, the photonic density 
of states in hyperbolic metamaterials is very higlP. It 
results in the enhancement of the radiated power for a 
dipole radiator embedded in such a mediurrpEEl. Usually, 
this enhancement of radiation is used for optical sensing 
and locating the fluorescent nanoobjects. The gain in the 
radiated power at the frequency of spontaneous emission 
by the fluorescent object implies the increase of the de¬ 
cay rate. This increase is called Purcell factor and may 


attain two orders of magnitude for emitters embedded 
into a hyperbolic metamateriaP®. Initially, the Purcell 
effedP 101 was treated as an increase of the spontaneous 
decay rate of an emitter located in a resonant cavity. Re¬ 
cently, this concept was generalized to the modification 
of the decay rate in presence of any environment different 
from free spac^H 

Uniaxial hyperbolic metamaterials operating in the op¬ 
tical range are usually implemented as a stack of bilayers. 
One layer of the unit cell is a plasmonic metal, another 
layer is a transparent material (dielectric or semiconduc¬ 
tor). The thickness of the plasmonic layer is smaller than 
the skin-depth that makes the unit cell penetrable for 
propagating waves. The total thickness of the unit cell 
is much smaller than the wavelength in the transparent 
material. Another known implementation of hyperbolic 
metamaterial is a regular optically dense array of par¬ 
allel metal wires called wire medium (see e.g. irpHZl) 
The homogenization models known for both kinds of hy¬ 
perbolic metamaterials (see irP^tSi) point out to the 
possibility of a sharp qualitative change of their disper¬ 
sion regime versus frequency. When the frequency varies 
one of two principal components of the permittivity ten¬ 
sor (ey or £j_) may keep its sign whereas the sign of the 
other component changes. Then the hyperboloid in the 
reciprocal space transforms into a closed surface, e.g. an 
ellipsoid. This transition is similar to the known Lifshitz 
topological transitions in metalJ^ (when the Fermi sur¬ 
faces change their topology similarly). Therefore, this 
jump of the dispersion properties was entitled “topolog¬ 
ical phase transition” in RefP^. Topological transition 
may additionally enhance the spatial spectrum of propa¬ 
gating eigenmodes of the medium and, consequently, the 
gain in the spontaneous decay rate/radiated power of an 
emitter located in this medium. In fact, the hyperboloid 
though infinitely extended implies the cut-off for the ax¬ 
ial component of the wavevector. This cut-off at a given 
frequency w is equal to fco£j_, where ko = oj/c is the free- 






2 


space wave number. The spatial harmonics whose axial 
wavevector is below cut-off do not propagate. 

In RefP-^, a constitutive bilayer of the hyperbolic meta¬ 
material is formed by titanium oxide and silver nanolay¬ 
ers. At the topological transition frequency, e± of this 
metamaterial tends to zero in the framework of the ho¬ 
mogenization model (neglecting optical losses). Then the 
distance between two branches of the hyperbola squeezes 
and the cut-off vanishes. The spatial spectrum at this 
frequency is not only infinitely extent, it is unbounded. 
This obviously grants the maximal possible density of 
photonic states. Therefore, the Purcell factor calculated 
in Ref ! 16 achieved two orders of magnitude even for a 
quantum dot located on the surface of the hyperbolic 
metamaterial. Of course, optical losses decrease this ef¬ 
fect, however, in work 1 ^it was experimentally shown (in 
this case for a stack of silver-germanium bilayers) that 
even the radiative Purcell factor - the gain in power ra¬ 
diated to free space one achieves one order of magnitude 
for emitters located on the surface. Further, in RefP^one 
presented a general theoretical study of the phenomenon 
of unbounded spatial spectrum where the possible mag¬ 
netic response (effective permeability tensor) was also in¬ 
volved. 

We think that the strong Purcell effect obtained at the 
topological transition of a stacked hyperbolic metamate¬ 
rial does not close the chapter. We can study the other 
realization of a hyperbolic metamaterial, wire medium 
which can be fabricated with free-standing part^Sl, alter¬ 
natively nanowires can partially stay in a liquid host"-®. 
This allows quantum emitters to easily emerge in the 
bulk. Moreover, wire medium has its own advantages 
in nanofabrication. Fraction ratio and parallelism of 
nanowires are important, but the strict periodicity is 
not crucial for material parameters of the wire medium, 
whereas the stacked implementation of hyperbolic meta¬ 
material is very demanding to the flatness and periodicity 
of constitutive nanolayers. 

In the available literature topological transitions in 
wire media, to our knowledge, has not been revealed, 
yet. This is not surprising because most popular wire me¬ 
dia are implemented of plasmonic (Ag or Au) nanowires. 
For plasmonic wire media with practical design parame¬ 
ters topological transitions, to our estimations, occur at 
frequencies where the homogenization model loses its va¬ 
lidity. Since the homogenization effective-medium model 
is analytical, it represents the best tool for searching the 
needed regime. Therefore, we inspected different materi¬ 
als with the use of the effective-medium model in order 
to find the needed regime. The finding was validated via 
the exact calculation of the Purcell factor. 

In this paper, we report the topological transition in a 
wire medium whose nanowires are prepared of LiTaC^. 
It occurs in the mid IR range, where the array period is 
much smaller than the wavelength and the homogeniza¬ 
tion model keeps adequate. Both hyperbolic and elliptic 
dispersion regimes are combined at the transition fre¬ 
quency, and the corresponding dispersion branches are 


connected that offers the unbounded spatial spectrum. 
Though optical losses in such the wire medium are sub¬ 
stantial they are not overcritical and the unbounded spa¬ 
tial spectrum manifests itself in the high Purcell factor. 
This “transitional” Purcell factor should be, evidently, 
higher than that obtained in Ref . 21 for a wire medium op¬ 
erating in the “standard” hyperbolic regime (with a cut¬ 
off). Since the spatial spectrum is unbounded over the 
axial wavevector the axial dipole should generate more 
radiation than the transverse one. This unusual expec¬ 
tation diverges from all previously known results. The 
Purcell factor of “standard” (hyperbolic) wire media at¬ 
tains two orders of magnitude for a tr ansverse dipole, 
whereas for an axial dipole it is smalleJ^ 7 * 21 * 22 ( In the 
transition regime the Purcell factor of a transverse dipole 
keeps the same as in the hyperbolic regime, whereas for 
an axial dipole it exceeds two orders of magnitude. 

Our paper is organized as follows: in Sec. [TTJ we intro¬ 
duce the approach how to achieve the unbounded spatial 
spectrum in a wire medium. Section |III| presents the 
numerical results for the Purcell factor and gives some 
discussions on their interpretation. Finally, Sec. [IV] con¬ 
cludes the paper. 


II. THEORY 

Electromagnetic properties of a general wire medium 
shown in Fig. |T] is described through a homogeneous ef¬ 
fective permittivity tensor as 
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FIG. 1. Long parallel wires which are set in a square lattice. 
b and a are the lattice constant and the radius of a dielectric 
rod, respectively, and £h represents the host medium. 
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the effective tensor are given bjfm 
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Here, £h and e r are the permittivity of the host medium 
and the wires, respectively, and f v denotes the volume 
fraction. In this paper, we assume that the host medium 
is air, i.e. £h = 1. In addition, ho is the free-space wave 
number, f3 is the axial component of the wavevector in the 
wire medium, and k p is called the plasma wave-number: 
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Since the wire medium is a uniaxial, two modes exist in 
it. The first mode is ordinary transverse-electric (TE) 
wave, and the second one is extraordinary transverse- 
magnetic (TM) wave. In this work, we are interested in 
the extraordinary wave because the hyperbolic dispersion 
corresponds only to this mode. Solving Maxwell’s equa¬ 
tions, the dispersion relation of extraordinary wave can 
be expressed a:P^ 




( 4 ) 


where q is the transverse component of the wavevector. 
Here, en is a function of /3. 

For simplicity, let us assume that the wires are lossless. 
If their volume fraction is small enough and the relative 
dielectric constant of the wire material is not very close 
to —1, the transverse component £j_ is positive and ap¬ 
proximately equals to unity (in the general case - to that 
of the host dielectric material). Equation[4]indicates that 
if £|| < 0, a hyperbolic dispersion arises, and if £m > 0, an 
elliptical dispersion arises. Since hyperbolic dispersion 
implies propagation of waves with high values of both 
axial and transverse wavevectors, the evanescent waves 
usually produced by a subwavelength dipole source are 
converted into propagating waves. This results in the 
high Purcell factor, as it was already mentioned in the 
introduction. 

As it is seen from Eq. [4j there is a cut-off for axial 
component of the wavevector: /3 > ko^/¥±- To increase 
the Purcell factor we need to eliminate this cut-off. For 
/3 < k 0y /e±, based on Eq. Bl the propagation becomes 
possible if £j_ > 0 (elliptical dispersion). We need to 
unify both elliptical and hyperbolic dispersions at the 
same frequency so that these branches touch each other 
at the transition point f) = koy/e±. It is possible because 
the axial component of permittivity in Eq. [4] depends on 
the wavevector. 


If we substitute Eq. [2]in Eq. 0 the transversal wavevec¬ 
tor can be expressed as 

q 2 = - [/3 2 - Ri] [/3 2 ~ R 2 } , 2 f f v(1 ~f;, ) 2 
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Ri = &0 £ -L- 


If the roots, i?i and i ?2 are equivalent and / v (l — £ r ) > 0, 
the value q 2 is always positive which results in propagat¬ 
ing mode for any /3. The cut-off of spatial spectrum is 
removed at this frequency, and we have 


k 2 

/ v (l - £ r ) = -TVT2 ■ 

k 2 + (£_l -l)^o 
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Since the plasma wave number is very large compared to 
y/ £j_ — 1*0 (fcp ^ (£_L ~ l)fco), the needed permittivity 
of the wire material can be approximated as 


£ r = 1 • (8) 

Jv 

Condition expressed in Eqs. JT] and [8] results in £y = 0 at 
the point /? = ko^/e± in the iso-frequency surface. Based 
on Eq. |2j it can be easily shown that for fi < koy/e±, £y 
is positive (corresponding to nearly elliptical dispersion), 
and for /3 > koy/e±, £m is smaller than zero (correspond¬ 
ing to nearly hyperbolic dispersion). This topological 
transition at the desired frequency removes the cut-off 
in the axial wavevector. Since the permittivity tensor 
depends on the wavevector, we do not obtain exactly hy- 
perboloidal and ellipsoidal parts of the dispersion surface. 
Therefore, instead of hyperboloid and ellipsoid, we prefer 
to use more exact terms: open surface and closed surface, 
respec tively . 

Fig. |2(a)| shows the axial component of the effective 
permittivity versus the normalized axial wavevector at 
the optimized wavelength in which the topological tran¬ 
sition happens. In this example the fraction of wires 
/ v = 0.05 that makes the relative dielectric constant 
of the wire material equal to £ r = —19 in according 
to Eq. [8] The figure has specified the closed-surface 
(positive-permittivity) and the open-surface (negative- 
permittivity) regimes. The transition happens at /36 = 
±1 where £||(/3) = 0. The corresponding iso-frequency 
curve is shown in Fig. |2(b)[ Different iso-frequency sur¬ 
faces for several values of relative dielectric constant of 
the wires are shown in Fig. [3j It is easily seen that the 
optimal regime corresponds to the blue solid curve rep¬ 
resenting £ r = —19 when there is no spatial-frequency 
cut-off. When £ r = —23, the closed surface does not ex¬ 
ist and the cut-off is present. Reducing gradually the 
absolute value of the relative dielectric constant causes 
the positive-permittivity regime to appear. Both parts 
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(b) 

FIG. 2. (a)-Axial component of effective permittivity with 
respect to normalized wavevector and (b)-corresponding iso¬ 
frequency contour (cross section of the dispersion surface). 
Calculations are related to the following values: f v = 0.05, 
and A/6 = 6.6 where A is the free-space wavelength. 


of the iso-frequency surface touch each other at the point 
£|| = 0 if e r = —19. 

So, the topological transition in a wire medium whose 
wires have relative permittivity close to —19 is possible in 
the frequency range where the non-local homogenization 
model® is applicable. This model has been recognized 
for its high accuracy 12 . At the corresponding frequency, 
the spatial spectrum is unbounded (no cut-off as for the 
stacked implementation of the hyperbolic metamaterial 
at its topological transitiorf®). However, this effect hap¬ 
pens not because the distance between two hyperbolic 
branches squeezes as in®. In our case both closed-surface 
and open-surface branches exist in the reciprocal space 
together and touch one another. 

Notice, that such the phenomenon is, in principle, 
known in the theory of metamaterials with spatial disper¬ 
sion. In work® one has suggested and theoretically stud¬ 
ied a spatially dispersive metamaterial for which both 
hyperboloid and ellipsoid surfaces exist in the recipro¬ 
cal space simultaneously (beyond a topological transi¬ 
tion) . The suggested metamaterial in® is an ultimately 
anisotropic strictly regular lattice of uniaxial strongly res¬ 
onant inclusions. In optics such an array could be imple¬ 
mented from very thin and substantially long Ag or Au 
nanorods operating at their plasmon resonance. Varying 
the geometric parameters one may engineer the regime of 



(3b 

(b) 


FIG. 3. Iso-frequency contours for different values of the rel¬ 
ative dielectric constant of the wires. Calculations are related 
to the following values: f v = 0.05, and A/6 = 6.6. The curves 
based on their colors correspond to red-e r = —23, black¬ 
er = —20.75, blue-e r = —19, green-e r = —17 and magneta- 
£r = -15.35. 


unbounded spatial spectrum when the open and closed 
branches of the dispersion surface would be touching or 
intersecting. 

However, the metamaterial® is challenging for practi¬ 
cal realization at optical frequencies, and its experimental 
verification in the near future is problematic. Meanwhile, 
the spatial dispersion which is needed to unify the elliptic 
and hyperbolic dispersion branches turns out to be possi¬ 
ble in a wire medium. In the next section, we consider an 
implementation of the wire medium of a polaritonic crys¬ 
talline material. Aligned nanowires of polaritonic materi¬ 
als grown in a solid matrix have been reported since 2009 
where the first paper on them was published^ 5 . Though 
for the instance polaritonic wire media with free-standing 
parts are not known, there are no basic technological lim¬ 
itations, and we hope that they will be reported soon. 
Moreover, the topological transition with needed spatial 
dispersion is not a prerogative of polaritonic crystals. We 
plan to find the needed regime for polycrystal media more 
popular among nanotechnologists such as SiC or GaN 
which possess the needed modest negative relative per¬ 
mittivity in the mid IR range. 
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III. RESULTS AND DISCUSSION 


and FPp = RPf- 


To confirm our expectation for the Purcell factor at 
the topological phase transition, we performed full-wave 
simulations employing the CST Microwave Studio sim¬ 
ulator. For reliability, some simulation results were also 
reproduced using the HFSS software. The sample of wire 
medium in our simulations is cubic, and the volume frac¬ 
tion of lithium tantalate is equal to / v = 0.0804. This 
value may correspond, for example to the wire radius 
a = 32 nm and period b = 200 nm. Based on Eq. [8j 
the relative dielectric constant of the wires should have 
Re(e r ) = —11.4 and |Im(e r )| -C |Re(e r )| that allows ey ~ 
0 at the point (3 = koy/e± and subsequently the tran¬ 
sition from negative-permittivity regime (Refey) < 0) 
to positive-permittivity regime (Re(ey) > 0). Because 
Re(e r ) is negative, the transition is in principle possi¬ 
ble for metal nanowires within visible range. However, 
this option requires either very tiny nanowires which 
will be challenging for fabrication or more sparse ar¬ 
ray for which the transition occurs at the wavelength 
beyond the domain of applicability of the homogeniza¬ 
tion model. Another option is using polaritonic ma¬ 
terials, e.g. lithium tantalate operating in the near- 
IR or rnid-IR ranges. The dielectric function of such 
materiales is provided by the Drude-Lorentz modeP^J 
e r = £oo[l + (w£ — w t)/( w t — w 2 + 7)] where for 

LiTaC >3 we have wt/27t = 26.7 THz, = 46.9 

THz, 7/277 = 0.94 THz, and £<*, = 13. 4 2 '®. The value 
e r = —11.4 corresponds to the permittivity of lithium 
tantalate at about / = 39 THz. 

The emitter is chosen to be a sub-wavelength electric 
dipole. To have studied different orientations and loca¬ 
tions of this dipole inside the wire medium sample. Here, 
we show the results for both longitudinal (dipole paral¬ 
lel to nanowires and optical axis) and transversal (dipole 
perpendicular to optical axis) orientations and compare 
the results. In both cases the dipoles are located sym¬ 
metrically in the gap between nanowires at the center of 
the wire medium sample, as it is shown in Fig. |4j 

We calculated two types of the Purcell factor. One 
is called full Purcell factor ( FPp) which shows the total 
enhancement of the radiated power of the dipole. It can 
be found through the real part of the input impedance 
of the lossless Hertzian dipole which is equal to its radi¬ 
ation resistance. Keeping the same current in the dipole 
in the presence of the metamaterial sample and in its ab¬ 
sence the ratio of radiation resistances delivers the gain 
in the radiated power FPp. The second type of the Pur¬ 
cell factor is called radiative Purcell factor {PPp). It 
determines the enhancement of the power radiated by 
the same dipole moment into free space. RPf is calcu¬ 
lated via the power flux integrated in the far zone in the 
presence of the sample and in its absence. An important 
check of the reliability of our simulations was done nulli¬ 
fying the optical losses in lithium tantalate “by hands”. 
Then the absorption of radiated power in the sample, i.e. 
the non-radiative part of the full Purcell factor vanishes 


A. Longitudinal dipole emitter 

The geometry and orientation of the dipole is shown 
in Fig. [4} 



FIG. 4. Geometry and orientation of the sub-wavelength elec¬ 
tric dipole source in the metamaterial sample. 


Figure 5(a) depicts the full Purcell factor versus fre¬ 


quency. To be sure that our effect is not distorted by 
dimensional resonances, we simulated the cubic sample 
with five sizes of the edge, corresponding to 10 x 10, 
12 x 12, 14 x 14, 16 x 16 and 18 x 18 nanowires, re¬ 
spectively (the internal structure of the wire medium 
was kept the same, that means the length of nanowires 
varying from 1.8 microns to 3.4 microns). As it is seen, 
the dimensional resonances do not enter the selected fre¬ 
quency range i.e. the full Purcell factor in all cases em¬ 
ulates the infinite wire medium. About the frequency 
f = 38 THz, the full Purcell factor has its maximum 
value, whereas the homogenization model predicts the 
maximum at / = 39 THz. This small error is not surpris¬ 
ing. First, the non-local homogenization model though 
most accurate of known effective-medium models of wire 
media is still an analytical approximation. Second, the 
whole theory of Sec. [H] neglects optical losses, though 
they are quite substantial. Notice, that these hi gh lo sses 
make the resonance of the Purcell factor in Fig. |5(a) not 
very pronounced. Indeed, in a lossy medium the topo¬ 
logical transition cannot be sharp. It happens not at a 
single frequency but across a certain frequency interval. 
Another reason why the resonance in Fig. |5(a)| is not suf¬ 
ficiently pronounced is the logarithmic scale of the plot. 

Whereas the full Purcell factor can exceed 600 at the 
transition frequency, the radiative Purcell factor shown in 
Fig. 5(b) | is not very high. However, it is still remarkable 
at the resonance. The global maximum was achieved 
for the sample of 10 x 10 nanowires, and the maximum 
occurs at nearly / = 37 THz. If the size of the sample is 
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FIG. 5. Purcell factor for different sizes of lithium tantalate 
cubic wire medium. (a)-Full Purcell factor. (b)-Radiative 
Purcell factor. Here a = 32 nm and b = 200 nm. 


FIG. 6. Purcell factor for different radii of wires in lithium 
tantalate cubic wire medium including 10 x 10 nanowires, (a)- 
Full Purcell factor. (b)-Radiative Purcell factor. Here b = 
200 nm. 


smaller, the sample becomes mesoscopic (the full Purcell 
factor starts to feel the sample size), and on the other 
hand, if the size of the sample is larger, the maximum 
of radiative Purcell factor decreases monotonously versus 
the size which means that the internal reflections have no 
visible impact. At the resonant frequency, the reflection 
at the interface with free space is not significant because 
the wave impedance of the sample is relatively matched 
to the free-space wave impedance. This matching allows 
the electromagnetic wave to exit from the sample. The 
smallness of the radiative Purcell factor compared to the 
full one can result from the attenuation of the radiation 
in the sample. 

As to other frequencies, around the topological phase 
transition range, the internal reflections become impor¬ 
tant due to the wave impedance mismatch. Close to 
the resonant frequency, the radiative Purcell factor drops 
dramatically because the metamaterial sample experi¬ 
ences the strong wave impedance mismatching at the 
interface. This mismatch definitely implies the stand¬ 
ing waves, however, the dimensional resonances do not 


arise due to the high attenuation. At low frequencies 
where the open-surface regime only exists, a local maxi¬ 
mum is seen that is most visible for the sample of 10 x 10 
nanowires. This local maximum at / = 32 THz is also 
due to the impedance matching since it exists for all sizes 
of the sample. 

Figures 6(a) and |6(b)| compare the full Purcell factor 
and radiative Purcell factor, respectively, for several radii 
of nanowires while the lattice constant is fixed. Accord¬ 
ing to Eq. [8j by changing the radius i.e. the volume 
fraction, we change e r and, consequently, shift the fre¬ 
quency of the topological transition. This frequency shift 
is clearly seen in the figures for both full and radiative 
Purcell factors. 


B. Transversal dipole emitter 

In this subsection, we report the results of similar cal¬ 
culations for the transverse dipole (with the same abso- 
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lute value of the dipole moment as above) and compare 
to the previous ones obtained for an axial one. For bet¬ 
ter clarity of our plots, we show the results only for two 
cubic samples - 10 x 10 and 12 x 12 nanowires. 



(a) 



Frequency [THz] 

(b) 


FIG. 7. Purcell factor for two different orientations of the 
dipole. Two samples of lithium tantalate cubic wire medium 
including 10 x 10 and 12 x 12 nanowires are considered, (a)- 
Full Purcell factor and (b)-radiative Purcell factor. Here, 
a = 32 nm and b = 200 nm. 


Figure 7(a) shows that the transverse dipole does not 
feel the topological transition. Definitely, the radiation 
decreases above / = 40 THz, where the dispersion sur¬ 
face is closed. However, at the frquency range where the 
topological transition happens for the axial dipole, the 
Purcell factor for the transverse one is the same as be¬ 
low - in the “standard” open-surface regime. This fully 
agrees with our theoretical expectations. The Purcell fac¬ 
tor of the transverse dipole turns out to be smaller than 
that of the axial dipole (this result keeps the same for 
other sizes of the sample) at the transition. Especially, 
the insensitivity of the transverse dipole to the topologi¬ 
cal transition is seen in the plot for the radiative Purcell 


factor depicted in Figure [7(b) | The longitudinal (axial) 
dipole is much more efficient for radiation to free space 
than the transversal one. For it, the radiative Purcell fac¬ 
tor does not exceed 2 and equals unity in the range of the 
topological transition and around. Definitely, this is the 
result of the strong wave impedance mismatch. Since 
the material is anisotropic, the regime ey « 0 mimics 
the e-near-zero regime only for the radiation of the axial 
dipole (with dominating axial polarization), whereas the 
radiation of the transverse dipole almost does not feel 
the zero of ey and turns out to be confined inside the 
sample. Therefore the contrast between RPp of the axial 
and transverse dipoles at about the transition frequency 
is 10-fold. 


IV. CONCLUSIONS 

In this work, we have theoretically revealed and studied 
topological phase transition in wire medium which may 
occur in the frequency range where the homogenization 
model is still valid. We utilized the spatial dispersion 
inherent to wire medium in order to engineer the un¬ 
bounded spatial spectrum of propagating eigenmodes at 
the transition frequency. The open and closed parts of 
the dispersion surface in this regime touch each other at a 
special point in the reciprocal space of the lattice. At this 
point the real part of the axial component of the effec¬ 
tive permittivity is zero. We found the practical design 
parameters for lithium tantalate nanowires which allow 
to implement this regime. 

Due to this effect (topological phase transition over the 
wavevector axis at one certain frequency), we obtained 
by exact numerical simulations a dramatic increase of 
the radiated power of a sub-wavelength electric dipole 
located at the center of the finite-size sample of our wire 
medium. The highest Purcell factor corresponds to the 
axially oriented dipole, that represents the contrast to 
previously known results where the high Purcell factors of 
wire media were reported only for the transverse dipoles. 
Whereas the radiation of the axial dipole is resonant with 
the maximum at the transition frequency, the radiation 
of the transverse dipole is rather stable in the frequency 
range of open-surface dispersion and smoothly decreases 
at higher frequencies where the dispersion surface be¬ 
comes closed. We also investigated the radiative Purcell 
factor. It attains one order of magnitude for an axial 
dipole. 

Finally, since lithium tantalate is a lossy material, and 
the radiative Purcell factor we have reported is modest, 
we plan to reveal the same regime for wire medium of 
SiC or GaN, where RPf may approach more closely to 
the resonant value of the full Purcell factor e.g. attain 
two orders of magnitude. Then we will organize an ex¬ 
perimental verification of our effect. 
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